Atomic Fermi-Bose Mixtures in Inhomogeneous and Random Lattices: 
From Fermi Glass to Quantum Spin Glass and Quantum Percolation 
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We investigate atomic Fermi-Bose mixtures in inhomogeneous and random optical lattices in the 
limit of strong atom-atom interactions. We derive the effective Hamiltonian describing the dynamics 
of the system and discuss its low temperature physics. We demonstrate possibility of controlling 
the interactions at local level in inhomogeneous but regular lattices. Such a control leads to the 
achievement of Fermi glass, quantum Fermi spin glass, and quantum percolation regimes involving 
bare and/or composite fermions in random lattices. 
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Fermi-Bose (FB) mixtures attract considerable inter- 
est in the physics of ultra-cold atomic and molecular 
gases, comparable with the interest in molecular Bose- 
Einstein condensation (BEC) 0, or Bardeen-Cooper- 
Schrieffer (BCS) transition |2( in ultra-cold Fermi mix- 
tures. The reason for interest in FB systems is threefold. 
First, these are very fundamental systems that have no 
direct analogues in condensed matter physics. Second, 
these systems can be very efficiently cooled using sympa- 
thetic cooling down to very low temperatures T (of order 
of tens nK) 0, U 0- ■ Finally, their physics is extremely 
rich and not yet fully understood. 

FB mixtures have been intensively studied in traps , 
but the experimental observation of the superfluid (SF) 
to Mott-insulator (MI) transition in bosonic gases |14j . 
predicted in Ref. [lfij . has triggered the interest in the 
physics of FB mixtures in optical lattices Under ap- 
propriate conditions such mixtures can be well described 
by the Fermi-Bose Hubbard model (FBH) . A partic- 
ularly appealing feature of the FBH model is the possibil- 
ity to produce novel quan tum phases , fermion-boson 
induced superfluidity [ijj, and composite fermions, which 
for attractive (repulsive) interactions between fermions 
and bosons, are formed by a fermion and bosons (bosonic 
holes) as first shown in {5(| (see also [2lt Eif ). 

FB mixtures in the limit of strong atom-atom interac- 
tions (strong coupling regime) show a very rich variety of 
quantum phases in periodic optical lattices |'2l| . They in- 
clude the mentioned composite fermions, and range from 
a normal Fermi liquid, a density wave, a superfluid liq- 
uid, to an insulator with fermionic domains. The phase 
diagram of the system has been determined in Ref. [2^] 
by means of mean- field theory jiij . These studies have 
been generalized recently to inhomogeneous lattices [2^] 
to include the effects of the lattice and of a possible trap 
potential. So far, only the case of strong interactions and 
vanishing hopping has been considered. 

In the present Letter we study the low temperature 
physics of FB mixtures in optical lattices with local and 



random inhomogeneities in the strong interactions limit 
but including tunneling as a perturbation. We show that 
interactions and tunneling may be controlled at the local 
level in inhomogeneous lattices [2fij . This control gives 
access to a wide variety of regimes and we derive the cor- 
responding effective Hamiltonians. We then show how 
to achieve Fermi glass, fermionic spin glass and quan- 
tum percolation regimes involving bare and/or composite 
fermions in random lattices. 

We consider a sample of ultra-cold bosonic and (polar- 
ized) fermionic atoms (e.g. 7 Li- 6 Li or 87 Rb- 40 K) trapped 
in an optical lattice. At low temperature, the atoms oc- 
cupy only the lowest energy band and it is convenient to 
work in the corresponding Wannier basis . Note that 
a fermion number N-p strictly smaller than the number 
of lattice sites N is required here. The Hamiltonian of 
the system reads [ljl |27| : 



#fbh = - X)( Jb6 J 6 J + Jf -^ + kc ') 
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where b\, bj, fj, fj are the bosonic and fermionic 
creation-annihilation operators, rii — b\bi and m t = f} fi- 
The FBH model describes: i) nearest neighbor (n.n.) bo- 
son (fermion) hopping, with an associated negative en- 
ergy, — Jb (— Jf); ii) on-site repulsive boson-boson in- 
teractions with an energy V; iii) on-site boson-fermion 
interactions with an energy U, which is positive (nega- 
tive) for repulsive (attractive) interactions, iv) and, fi- 
nally, interactions with the optical potential, with ener- 
gies nf and /if. In the following, we shall consider only 
the case Jb = Jf — J and the regime of strong interac- 
tions, namely V,U ^> J. 

In a periodic optical lattice, fif ' F is simply the (bosons 
or fermions) chemical potential and is independent of the 
site i. It is, however, possible to add a laser field inde- 
pendent of the lattice to modify the depths of the optical 
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potential wells in a site-dependent way |28j . In this case, 
the local potential depth has to be added to fxf ,F , which 
may now be inhomogeneous. If the added field is periodic 
and if the spatial period is commensurate with the lattice 
period, fa ' F is periodic; if the spatial periods are incom- 
mensurate, fa ' is quasi-periodic. One can also add a 
random speckle field, so that [i i ' is random. Experi- 
mental techniques offer full possibilities to control such 
periodic, quasi-periodic, or disordered fa' F |29| . Note, 
that the additional inhomogeneous potential might, but 
does not have to, act equally on both atomic species. 
Here, we will study the case /if = 0, fa 3 = faV. 

In Ref. [HI we have used the method of degenerate 
second order perturbation theory to derive an effective 
Hamiltonian by projecting the wave function onto the 
multiply degenerate ground state of the system in the ab- 
sence of tunneling. This can be extended to the present 
situation, where there are very many states with similar 
energies. It is thus reasonable to project the wave func- 
tion on the manifold of 'ground states'. These states are 
local minima of energy, since at least some of hopping 
acts increase their energy by V or \U\. 

Let us consider first J = 0, and the case < fa < 1. 
In the absence of a fermion one expects one boson per 
site, i.e. n,- = 1 j2lj. We shall consider here only the case 
of repulsive interactions, i.e. a = U/V > 0. 

It is useful to divide the sites into : i) ^4-sites, for which 
fa — a > 0, and fermions do not push bosons out, and 
ii) -B-sites, for which fa — a < 0, and the fermion pushes 
the boson out forming a composite fcrmion-bosonic hole. 
Energetically, the second situation is favorable, so for a 
given set of N A of A-, and N B = N - N A of B-sites, 
the fermions at very low temperature will first occupy 
the B-sites until N-p = N B , and then they will start 
to occupy the A-sites. We construct the corresponding 
projector operators P, Q = 1 — P, which depend on 
N A and Np. The operator P describes the projection 
onto the manifold of quasi-degenerated states in which 
the fermions occupy the P-sites stripped of bosons and 
some ^4-sites only if iVp > N B . In this case there is a bo- 
son in any A-site and if Np < N B there is also a boson 
in the P-sites which do not contain fermions. We use 
the formalism of second order time dependent perturba- 
tion theory |32| . and project the Schrodinger equation, 
ihdt\ip(t)) = (Ha + Hi)\ip(t)), onto the manifold of states 
spanning P. The "zeroth-order" part Hq contains the 
atomic interactions and terms proportional to the chem- 
ical potential and commutes with P. H\ represents the 
tunneling terms. The effective equation for \ipp) = P\ip) 
reads then ihd t \ipp{t)) = H e ff\ijjp(t)} where 

(out\H cff \in) = (out\H \in) + (out\P HxP\in) (2) 
-\{out\PHiQ ( 1 + 1 ) QPiPM. 



The effective Hamiltonian P e // has the form 

Pes = ^[-(dijF^Fj + h.c.) + KijMiMj] + ^ ftAfi, 
(y) * 

(3) 

where Pi, F- are the corresponding (composite) fermionic 
annihilation and creation operators, and Mi = FjFi. 
The hopping amplitudes dij and the n.n. couplings 
(which might be repulsive (> 0) or attractive (< 0)) are 
of order of J 2 /V. The values of the couplings depend 
on a, fli, flj, and J, and have to be determined care- 
fully for different cases, as discussed below. Note, how- 
ever, that the hopping i — > j, or back causes the energy 
change ±(Ay = /i; — jij) in units of V, i.e is highly non- 
resonant and inefficient for Ay ~ 1; it will first lead to 
jump rates of order 0(J 4 /V 3 ). Additionally, composite 
fermions may feel the local energy jli. 

I) All sites are of type B. In this case we have a gas of 
composites flowing within the MI with 1 boson per site. 
The couplings are 




The chemical potential fa/V ~ fa up to corrections of 
order 0(J 2 /V). The hopping amplitudes dij are for this 
case always positive, although may vary quite signifi- 
cantly with disorder, especially when Ay ~ a. As shown 
in Fig. 1, for a > 1, Kij < and we deal with attractive 
(although random) interactions. For a < 1, but close 
to 1, Kij might take positive or negative values for Ay- 
small or Ay ~ a. In this case the qualitative character 
of interactions is controlled by inhomogeneity. 
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Figure 1: Nearest neighbor couplings Kij as a function of 
Ay. Solid line: Coupling in case I, with ai = 0.93. Dashed 
line: Same expression with ai = 1.07. Dashed dotted line: 
Coupling in case II with an = 0.03. 

At low temperature the physics of the system depends 
on the relation between fa's and a. For small inhomo- 
geneities, we may neglect the contributions of Ay to 
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dij ~ d and K^ ~ if, and keep only the leading disor- 
der contribution in /i,-. Note, that the latter contribution 
is relevant in ID and 2D leading to Anderson localiza- 
tion of single particles [33j]. When K <C d the system 
will then be in the Fermi glass phase, i.e. Anderson lo- 
calized (and many-body corrected) single particle states 
will be occupied according to the Fermi-Dirac rules [34| . 
For repulsive interactions and K d, the ground state 
will be a Mott insulator for large enough filling factors. 
In particular, for filling factor 1/2 one expects to form 
a checker-board phase. For intermediate values of K/d 
delocalized metallic phases with enhanced persistent cur- 
rents are possible |3Jjj. Similarly, for attractive interac- 
tions (K < 0) and \K\ < d one expects competition 
between pairing of fermions and disorder. In this case 
for \K\ d, the fermions will form a domain insulator. 

Another interesting limit is when |Ay| ~ a ~ 1. 
The tunneling becomes then non-resonant and negligi- 
ble, while the couplings fluctuate strongly. We end 
up then with the (fermionic) Ising spin glass model 
described by the Edwards- Anderson model |37j : 



E-A 



= 1^2 K ij s i s j + ^2 Mi^iA 



(0) 



with Si = 2Mi — 1 = ±1. The above Hamiltonian is 
well approximated by a random one with Gaussian and 
independent distributions for and fli/2 with mean 

(H), and variances K (h), respectively. In this limit 
the system may be used to study various open questions 
of spin glass physics, concerning the nature of ordering 
(Parisi's ^ versus "droplet" picture [3jJ |3JJ), broken 
symmetry and dynamics in classical (in absence of hop- 
ping) and quantum (with small, but nevertheless present 
hopping) spin glasses |2J,|4C|. The predictions of Parisi's 
mean field theory for the model © can be obtained by 
replacing the model by the corresponding Sherrington- 
Kirkpatrick model, and employing the standard method 
of replica trick [3jj . The calculations differ from the stan- 
dard ones in that the constraint of fixed mean number 
of fermions is applied, and one deals simultaneously with 
random couplings and "magnetic fields" /2j. Following 
de Almeida and Thouless (A-T) approach j^H, we ob- 
tain the A-T surface separating the stable paramagnetic 
state from the "true" spin glass state, characterized by 
replica symmetry breaking, and ultrametrically arranged 
ground states. The paramagnetic state is stable for 
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where q = ( ( tanh 2 



xy/K 2 q+h 2 +H 



(7) 



, the con- 



straint is in = yytanh ^ +H jjj ' w ^ m ~ 

2Ny/N — 1 and {{.)) x denotes averaging over normally 



distributed random variable x which represents disorder 
within the replica method . Note, that according to 
the predictions of the alternative "droplet" model |3S| - 
applied to ©, no A-T surface is expected to exist. 

II) All sites are of type A. In this case a < 1, and we 
have a gas of bare fermions flowing over the MI with 1 
boson per site. The coefficients are 



J, Ki. 



J2_ 

V 



(8) 
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and \ii ~ up to corrections of order 0(J 2 /V). The 
couplings Kij are positive, and for a ~ 0, Kij ~ 0{a 2 ), 
and both the repulsive interactions, and disorder are very 
weak, leading to a Fermi liquid behavior at low T. For 
finite a, and Ajj ~ 1 — a, however, the fluctuations of 
K^ might be quite large. Note, that for a ~ 1, this will 
occur even for small disorder. Assuming for simplicity 
that K^ take either very large, or zero value, we see that 
the physics of bond percolation W1 U pl a Y a r °l e - The 
bonds will form a "weak" and "strong" clusters, each of 
which may be percolating. The fermions will hope freely 
in the "weak" cluster; only one fermion per bond will be 
allowed in the "strong" cluster. 

Ill) Both N A and N B of order N/2. In this case 
the physics of site percolation ^3] wm be relevant. If 
N-p < N B the composite fermions will move within a 
cluster of B sites. When N B is above the classical per- 
colation threshold, this cluster will be percolating. The 
expressions Eq. @ and Eq. (J5J will still be valid, except 
that they will connect only the i?-sites. 

The physics of the system will be similar as in the case 
I), but it will occur now on the percolating cluster. For 
small disorder, and K <C d the system will be in a Fermi 
glass phase in which the interplay between the Ander- 
son localization of single particles due to fluctuations of 
Hi and quantum percolation effects, that is randomness 
of the B-sites cluster, will occur. For repulsive inter- 
actions and K ^> d, the ground state will be a Mott 
insulator on the cluster for large filling factors. It is an 
open question whether the delocalized metallic phases 
with enhanced persistent current of the kind discussed in 
Ref. (3^1 might exist in this case. Similarly, it is an open 
question whether for attractive interactions (K < 0) and 
\K\ < d pairing of (perhaps localized) fermions will take 
place. If | if | ^> d, we expect the fermions to form a 
domain insulator on the cluster. 

In the "spin-glass" limit Ay ~ a ~ 1, we deal with 
the Edwards- Anderson spin glass on the cluster. Such 
systems are of interest in condensed matter physics |43j . 
and again questions connected to the nature of spin glass 
ordering may be studied in this case. 

When N F > A B , all B-sites will be filled, and the 
physics will occur on the cluster of A sites. For a ~ 0, 
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we shall deal with a gas with very weak repulsive interac- 
tions, and no significant disorder on the random cluster. 
This is an ideal test ground to study quantum percola- 
tion at low T. For finite a, and Ajj ~ 1 — a, the interplay 
between the fluctuating repulsive Kij 's and quantum per- 
colation might be studied. 

Summarizing, we have studied atomic Fermi-Bose mix- 
tures in optical lattices in the strong interaction limit, 
and in the presence of an inhomogeneous, or random on- 
site potential. We have derived the effective Hamiltonian 
describing the low temperature physics of the system, 
and shown that an inhomogeneous potential may be effi- 
ciently used to control the nature and strength of (boson 
mediated) interactions in the system. Using a random 
potential, one is able to control the system in such a 
way that its physics corresponds to a whole variety of 
quantum disordered systems: Fermi glass, fermionic spin 
glass, and quantum percolation systems. 
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